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itCM^-  ABSTRACT 

Finite- amplitude  standing  wave  effects  in  air  at  ambient 
conditions  contained  in  rigid-walled  cylindrical  tubes  with  large 
length  to  diameter  ratios  were  experimentally  investigated.  These 
results  were  compared  to  a  perturbation  solution  of  the  one-dimensional 
non- linear  acoustic  wave  equation  which  incorporates  the  dissipative 
effects  due  to  viscous  and  thermal  energy  losses  at  the  walls.   The 
lowest  resonance  frequencies  of  the  tubes  ranged  from  62.5  Hz  to 
1  kHz,  and  sound  pressure  levels  (based  on  the  fundamental)  ranged 
up  to  160  dB.   The  finite-amplitude  distortion  was  in  excellent 
agreement  almost  up  to  the  onset  of  shock.   A  detailed  investigation 
of  small  amplitude  attenuation  in  standing  wave  tubes  was  conducted 
and  compared  with  the  Kirchhoff  equations.   Agreement  within  one 
percent  was  obtained  when  consideration  was  given  to  the  numerical 
analysis  correction  and  end-effect  losses. 
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1.   Introduction 

It  has  been  known  for  a  long  time  that  for  the  case  of  a  plane 
elastic  wave  of  any  amplitude  in  a  nondissipative  medium  the  exact 
solution  to  the  equations  of  hydrodynamics  predicts  a  change  in  form 
as  the  wave  travels  (1,2,3,4,5,6).   In  going  one  step  further  and 
considering  that  the  medium  is  dissipative,  it  can  be  shown  that  only 
waves  of  relatively  large  amplitude  evidence  any  appreciable  change 
in  form  (4,7).   This  provides  a  good  definition  for  finite-amplitude 
waves.   Stated  another  way,  a  wave  would  be  of  finite-amplitude  when 
the  inclusion  of  second- order  terms  are  required  for  the  theory  to 
predict  a  change  in  shape  as  the  wave  propagates. 

An  important  contribution  beyond  the  work  done  by  the  early 
investigators  of  the  nineteenth  century  was  by  Fay  (8).   He  investigated 
the  change  of  type  in  a  periodic  finite-amplitude  plane  wave  of  infinite 
extent  in  a  viscous  medium.   He  formed  a  Fourier-series  solution  valid 
far  from  the  source  that  was  an  attenuated  saw-tooth  wave. 

In  general  the  more  recent  theoretical  approaches  to  finite- 
amplitude  effects  have  been  of  two  types.  Approximation  methods  using 
some  assumption  on  the  manner  of  harmonic  generation  and  absorption  in 
a  medium  have  been  used  (9,10,11,12^  and  approximate  solutions  to  the 
nonlinear  differential  equations  by  the  perturbation  technique  or  other 
means  have  been  investigated  (8,13,14,15,16). 

An  expression  for  the  spatial  attenuation  of  an  plane  saw-tooth 
wave  finite  in  extent  was  developed  by  Rudnick  (17,18)  using  a  "first- 
order"  theory.   Since  the  wave  was  confined,  consideration  was  given 
to  the  attenuation  due  to  the  walls  by  using  the  results  of  Kirchhoff 
and  Helmholtz. 
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Goldberg  (13,14)  undertook  a  solution  of  the  problem  of  plane 
waves  in  an  unbounded  medium  by  a  direct  perturbation  approach  on 
the  equations  of  motion.   He  obtained  a  number  of  second-order 
solutions  and  was  able  to  provide  conclusions  about  the  relative 
importance  of  nonlinearity  and  dissipation. 

The  perturbation  analysis  has  been  extended  to  the  sixth  order 
by  Keck  and  Beyer  (15)  for  periodic  plane  progressive  waves  of  infinite 
extent  in  a  viscous  medium.   The  wave  equation  is  the  same  one  treated 
by  Fay,  but  the  solution  is  valid  near  the  source. 

A  numerical  solution  for  approximately  the  same  problem  has  been 
obtained  by  Cook  (20).   An  exact  solution  in  parametric  form  was 
obtained  by  Blackstock  (21)  for  the  problem  of  a  finite-amplitude 
sound  wave  produced  by  arbitrary  motion  of  a  piston  in  a  lossless 
tube. 

Most  investigators,  like  the  ones  already  discussed,  have  confined 
their  studies  to  traveling  waves.   However,  a  few  have  given  consid- 
eration to  standing  waves. 

An  exact  solution  of  the  nonlinear  wave  equation  for  the  case  of 
finite-amplitude  standing  waves  in  a  closed  tube  as  well  as  for  one 
having  a  vibrating  piston  at  one  end  was  obtained  by  Keller  (22). 
Upon  assigning  a  value  of  -1  to  the  adiabatic  exponent^,  the  nonlinear 
wave  equation  becomes  linear,  and  this  fact  was  used  by  Keller  in  his 
solution.   These  results  "blow-up"  at  resonance  because  the  dissipative 
forces  are  neglected  in  the  solution. 

Betchov  (23)  constructs  a  solution  at  resonance  based  on  discontinuous 
linearized  solutions  and  a  secular  equation.   This  solution  is  interesting 
in  that  for  an  invisid  fluid  the  amplitude  at  resonance  is  found  to  be 
finite  and  to  be  determined  by  non-linear  effects.   Betchov  also  discusses 
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the  modifying  effect  of  wall  friction  on  his  solution  and  assumes 
that  this  effect  is  equivalent  to  a  body  force  proportional  to 
the  velocity. 

Saenger  and  Hudson  (24)  present  a  simple  analytical  description 
of  the  behavior  of  periodic  shocks  at  resonance  using  a  two  part 
solution  of  the  nonlinear  one  "dimensional  wave  equation.   A  point 
of  importance  in  their  solution  is  that  the  amplitude  of  the 
oscillations  remains  finite  only  because  of  the  action  of  shear 
viscosity  and  heat  conduction. 

Chester  (25)  has  predicted  asymmetries  in  finite  amplitude 
standing  waves  with  an  elaborate  theoretical  investigation  based 
on  the  nonlinear  acoustic  equations.   The  solution  includes  the 
effects  of  compressive  viscosity,  and  of  shear  viscosity  in  the 
boundary  layer  at  the  walls  of  the  tube;  however,  his  solutions 
are  not  amenable  to  detailed  experimental  comparison. 

Coppens  and  Sanders  (26)  developed  an  extension  of  Keek-Beyer 
perturbation  approach  that  included  both  the  bulk  losses  and  the 
wall  losses.   Their  results  provide  information  on  the  amplitudes 
and  phases  of  the  Fourier  components  of  the  distorted  pressure 
waveform.   The  theory  is  applied  to  the  problem  of  finite-amplitude 
standing  waves  in  a  cylindrical  resonance  tube  that  is  excited  by 
a  piston  at  one  end.   This  perturbation  approach  was  extended  by 
Ruff  (27), and  computer  programs  were  written  to  extract  more 
information  from  the  theory. 

The  purpose  of  this  research  was  to  conduct  a  detailed 
experimental  investigation  of  finite-amplitude  standing  waves  in 
air  contained  in  rigid-walled  tubes  to  verify  the  soundness  of  the 
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Coppens-Sanders  theory  as  extended  by  Ruff.   In  addition,  a  detailed 
investigation  of  the  attenuation  of  sound  in  rigid-walled  standing 
wave  tubes  with  attention  given  to  end-effect  losses  was  conducted 
to  check  the  validity  of  applying  the  Kirchhoff  attenuation  equations 
to  standing  waves. 
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2.   Theory 

2.1   General 

A  plane  elastic  wave  traveling  in  a  fluid  exhibits  a  change 
in  form  as  it  travels.   This  can  be  seen  from  the  nonlinear  properties 
of  the  equations  of  state  and  motion  which  follow.   This  theoretical 
development  is  extracted  from  that  of  Coppens  and  Sanders  and  is 
presented  here  to  aid  in  interpreting  the  experimental  work. 

The  adiabatic  equation  of  state  of  a  fluid 

can  be  represented  by  a  truncated  Taylor's  expansion  in  the  form 

f*(/f/pj(f-f>c)+     ( &/2poZ)  (jO-fif  (2) 

where 

p   =  acoustic  pressure 

O      -   instantaneous  density  at  a.  point 
/ 

P0     =  equilibrium  density 

B/A  =  parameter  of  non linearity. 
If  the  fluid  is  an  ideal  gas,  then  B/A  =  ~fr  - 1  where  0   is  the  ratio 
of  principal  specific  heats.   The  coefficient  A  can  be  expressed  as 

2. 


A  9fic 


*  ~<>  (3) 

where 

c   =  infinitesimal  amplitude  free  field  phase  velocity. 
The  equation  of  continuity  for  one  dimensional  wave  propagation 
can  be  written  in  Lagrangian  coordinates  as 
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where  £    is  the  particle  displacement  from  the  initial  position 
at  a. 

The  force  equation  can  be  written  in  a  slightly  generalized 
form  as 

/&  ^4£  -  ~$i<fi4'-Z*  (5) 

when  u  =  c^  >  is  the  particle  velocity  and  the  operator  FT  generates 
those  body  forces  in  Lagrangian  form  that  are  active  in  the  system 
of  interest. 

The  combination  of  Eqs .  2-5  yields 

3fe  -  tf"*§)*[i  -CM^O^f]  aft  =/>"#«.         (6) 

which  for  an  ideal  gas  becomes 

3t£  -q?0**4) _>/  3J  » /*"j 'a  .  (7) 

For  acoustic  processes  where  dj£"«/     Eq.  7  can  be  approximated  by 
Differentiation  of  Eq.  8  with  respect  to  time  yields 

(tf  +  oQJu  tf&t&S?  (9) 

where 

W  (10) 

Equations  6-9  have  been  the  basis  of  a  number  of  theoretical 
investigations  utilizing  successive  approximations  or  perturbation 
expansions  to  describe  finite-amplitude  effects  for  traveling  waves 
in  fluids  (with  or  without  viscous  dissipation)  (8,9,10,15,21,28). 
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Under  the  assumption  of  weak  dissipative  processes  and  with  the 

substitution  of  the  perturbation  series   s"  —  )    /      Eq.  9 

£         -*,  n  -V..  > 

becomes 

T.^at+A)^-  <k  (id 

where  the  sum  over  j  is  understood  to  encompass  just  those  frequencies 
appearing  in  fh  with 


0;,,   -  0  X>7=/  i  (lib) 

To  predict  the  finite-amplitude  effects  for  one-dimensional  wave 

propagation  in  a  fluid  system,  it  is  only  necessary  to  obtain  the 

linearized,  absorptive  wave  equation  for  the  system,  substitute  the 

resulting  expression  for  D   into  Eq.  11,  and  solve  for  the  frequency 

components . 

2.2   Rigid-Walled  Tubes 

Weston's  (29)  approach  to  the  propagation  of  sound  waves  in  gases 

contained  in  wide  tubes  can  be  used  in  the  formulation  of  D  .   For  the 

u 

wide  tube  the  sound  energy  is  diffused  evenly  over  the  tube  cross 
section,  and  the  adiabatic  motion  approximates  to  that  in  an  unbounded 
gas.   There  is,  however,  a  narrow  boundary  layer  at  the  wall,  analogous 
to  the  Prandtl  boundary  layer  in  continuous  flow,  where  viscous  and 
heat  conduction  processes  lead  to  an  energy  loss.   The  losses  in  the 
boundary  layer  at  the  wall  are  much  greater  than  the  bulk  losses. 
Weston's  equation  in  Eulerian  form  can  be  written  as 

where      -,  ,  _,  .  N  -    .-^   ,, 

f**(G/S)CZ<A^    -  1  (13) 
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^i  -    -yi  -[-(i^t)  C^'/a)  t      =  Kirchhoff  's  constant      (14) 

S  =  (uniform)  cross-sectional  area 
G  =  perimeter  of  area  S 

CO  =    (angular)  frequency  of  the  propagating  waveform 
V  =  kinematic  shear  viscosity  of  the  fluid 
y  =  thermal  diffusivity  or  thermometric  conductivity 

i  =  v^i: 


The  real  part  of  u  represents  the  mainstream  particle  velocity 

in  the  tube.   For  Eq.  12  to  be  valid  ^<T<5 .where ^fl is  the  boundary 

layer  thickness  given  approximately  by  5(~\>/oJ)2    (29).   Equation  12 
can  be  written  as 

tgu  +  CS/a>  )  4a4  a  -Sd?U  =0  (is) 

with   the  assumption  of  one -dimensional   flow. 
Then 

<$£  *  ( fy  /% )  4'24  "  4'  <V"  < l6  > 

where  ,j  .  represents  the  value  of  the  dissipation  parameter  ^  -  f  or  a 

frequency  qj.  . 

The  approximate  wave  equation  is  now  given  by 

E/°**4&    *^  )*^4&S£ll<#?*       (17) 

For  a  resonance  tube  of  uniform  cylindrical  cross-section  excited  at 

a  =  0  by  a  piston  with  the  same  cross  sectional  area  driven  with  constant 

acceleration  at  an  angular  frequency  CO     and  terminated  at  a  =  L  by  a 

plane  rigid  boundary, the  appropriate  boundary  conditions  can  be  written 

as 

Clu  016)  ~  #</«>)  smC&I-Q))  (is) 
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^lt(^t)^^Cqe)^^f'Xt)  ~o     m>/e    (19) 

The  phase  angle  (J   will  be  specified  later.   From  these  boundary 
conditions  Coppens  and  Sanders  present  a  solution  that  describes 
the  finite-amplitude  waveform.   Certain  quantities  developed  by 
Coppens  and  Sanders  are  important  in  any  experimental  investigation 
and  are  therefore  reproduced  here.   The  propagation  constant,  k,  and 
the  small  amplitude  attenuation  constant,  OC ,    are  given  by 

J,  »U)/o  *  (W/GK/rtS,)  (20) 


and 


(21) 


CjWill  be  discussed  in  greater  detail  in  Section  3.3. 
Various  frequency  relationships  in  terms  of  the  infinitesimal- 
amplitude  resonance  frequency,^  ,  are  given  by 

Ao5  =  CO  —COp  (22) 

for  the  frequency  differential  from  resonance  and 

C04    "&>;>    (!+£$,)  (23) 

for  the  frequency  at  which  each  harmonic  making  up  the  pressure  wave- 
form has  decreased  in  value  to  l/""\/2  of  its  maximum  value.   The 
dissipation  parameter  for  higher  order,  J,,,  ,  is  given  by 

fM  -  S,  /  y>*  (24) 

where  the  integer  n  represents  the  harmonic.   It  is  also  convenient 
to  introduce  the  peak  Mach  number  of  the  first-order  solution 

=  \Jif  /Co  (25) 
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3.     Experii     atal  Cons]  lerations 

Insight   into   the    finite-amplitude  case   can  be  gair        by  fir 
relating  the  the    infinitesimal-ampl]  -'i  Lew 

toward   isol    ting    the    experimentally  determinable   pai  „ 

the  expressi* 

U,      (Aa/a>)(sw/>*cCL  Hrt'Af4  (26) 

and 

f      '  <*$,Uf    c  9  (21) 

the  press?  :gid  end  of  the  tube  can  be  found  to  be 

•     r       qCo  4  fo)($ifih  ccL  Lsml&L)  $mcotm      (28) 

When  the  ace  I    : i-  a   of  the  piston  is  written  as 

j\r-:  t  )    a  4€   I  W'd0*-3)j  (29) 

it   can  be    setii    tl  piston   leads    the    first-order   pressure     at    the 

o 
rigid  end  of  tl       tubs     by   90     -    @1  .      The   first-order   p;  amplitude 


at   the  ria:      -  -  ttenuation   is   small 


!/      -   (^  ■    '<*>}{{ttLfi     ,  (30) 


and    the   pi  0,    becomes 

r  &  1 


..  ,,  ,    i/,)/CotL) 


■     .  ■   -<■  «."  <■  -  >-  /   .  31) 

At  resonance  Eqs.  30  and  31  become 

s/\Q/cCrm%  (32) 


and 


(33) 
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where  the  integer  m  represents  the  normal  mode  most  strongly  excited  by 
the  input  frequency.   When  finite-amplitude  effects  are  taken  into 
consideration,  the  pressure  is  modified  by  the  presence  of  the  higher 
harmonics  and  can  be  written  as 

It  is  convenient  to  introduce  the  resonance  parameter  Q  and  the  phase 

n 

parameter  9  defined  by 

am  -i  (**%>)/£  (35) 

and 

&*>  ®m    =(4~4n~  2^/^r  )/V; ,  d  (36) 

The  dependence  on  the  small  amplitude  attenuation  coefficient  ^Vis  best 
introduced  at  this  point  by  writing  the  dissipation  parameter  J/)J  in 
the  form 

£,  =  Z^^/cO/lt*  (37) 

through  the  use  of  Eqs.  24  and  21. 

If  0  and  0  are  rewritten  in  terms  of  <Xj  under  the  assumption  that 
the  bulk  losses  are  less  than  the  boundary  layer  losses,  then  the 
approximate  forms  are  given  by 

qma.^^c^%%/4^o>  (38) 


and 


\rift   %  ^    !    ->:      ~      r>&Jc0,    *CC0)t  (39) 


For  a  resonance  tube  with  a  diameter  of  1.125  inches,  a  length  of  6 
feet,  and  a  fundamental  resonance  frequency  of  approximately  100  Hz, 
the  bulk  losses  do  not  equal  the  boundary  layer  losses  until  the  input 
frequency  i»s  of  the  order  of  10  Hz.   This  would  be  in  the  range  of 
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4 
the  10  harmonic,  and  it  is  readily  seen  that  the  assumption  that  the 

boundary  layer  losses  dominate  is  valid. 

Equations  26-39  can  be  readily  used  for  a  comparison  with  the 

theory.   The  factors  A    Atand  c   are  known  for  the  system  to  be 

used,  and  the  quantities  A  and  COCtl/jJave   known  input  variables. 

0"}P0.3no <?0    can  be  determined  experimentally.   The  pressure  amplitude 

P   can  be  measured  directly,  and  0  can  be  determined  from  the 
nn  J  n 

frequency  dependence  of  the  pressure  measurements. 

Examination  of  certain  of  these  terms  reveals  some  important 

properties.   Q   is  maximum  above  <X*L,fcand  the  higher  the  subscript  n 

n  ~~ 

the  higher  the  maximizing  frequency.   With  constant  acceleration,  the 

higher  harmonics  have  their  maximum  amplitudes  at  successively  higher 

frequencies  with  narrower  peaking,  and  the  frequency  region  for  most 

favorable  finite-amplitude  excitation  lies  between  T/Jl  and  CO.  ^U)m(c/c)- 

The  angle  9  describes  the  difference  in  phase  between  the  piston 
n 

acceleration  at  a  =  o  and  the  pressure  at  a  =  L. 
3.1.   System  and  Equipment  Description 

A  diagram  of  the  system  used  in  this  experimental  investigation 
is  shown  in  Fig.  1.   The  resonance  tubes  contained  air  at  ambient 
temperature  and  pressure.   In  order  to  insure  as  nearly  as  possible 
ideal  boundary  conditions  and  to  reduce  the  effects  of  mechanical 
coupling  with  the  exciter,  the  resonance  tubes  were  fabricated  of 
steel  with  an  inside  diameter  of  2.250  inches  and  a  wall  thickness 
of  1.125  inches.   The  inside  surface  of  the  tube  wall  was  cleaned, 
but  no  other  attempt  was  made  to  improve  the  finish.   The  rigid  end 
was  terminated  with  an  aluminum  cap  2.0  inches  thick.   To  insure  no 
leaks,  the  surface  of  the  cap  and  the  end  of  the  tube  was  machined 
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and  polished  and  high  vacuum  silicone  grease  was  used  as  a  sealer. 
The  cap  was  held  in  place  by  five  studs  screwed  into  the  tube  wall. 
A  total  of  five  tubes  (identical  except  for  length)  were  used  to 
give  a  spread  in  fundamental  resonant  frequencies  in  excess  of  an 
order  of  magnitude.   The  lengths  of  the  tubes  were  7,  14,  and  28 
inches  and  6  and  9  feet.   The  six  foot  tube  was  used  for  the  major 
portion  of  the  investigation  and  for  comparison  with  the  computer 
predicted  results. 

The  piston  design  was  critical  for  several  reasons.   It  had  to 
be  heavy  enough  that  a  minimum  of  distortion  was  introduced  by  the 
force  exerted  upon  it  by  the  distorted  standing  wave,  yet  light  enough 
that  desired  acceleration  amplitudes  could  be.  obtained.   A  piston 
weight  of  approximately  0.8  pounds  was  used.   A  clearance  of 
approximately  0.006  inch  was  provided  between  the  seat  and  the  piston, 
and  an  air-tight  seal  was  obtained  by  the  use  of  two  0-rings  and 
silicone  grease.   The  0-ring  grooves  in  the  piston  were  designed  so 
that  the  neoprene  0-ring  would  be  under  compression  and  were  spaced 
to  provide  a  reservoir  for  the  grease.   A  smooth  seat  for  the  piston 
was  provided  by  machining  and  grinding  the  inside  of  the  tube.   The 
plane  surface  of  the  piston  was  machined  and  polished. 

The  M-B  Electronics  Model  EA  1500  exciter  due  to  its  versitility, 
small  size,  and  maximum  acceleration  rating  of  124  g's  was  chosen 
as  the  acoustic  driver.   For  the  production  of  intense  standing  waves 
in  a  resonance  tube,  large  displacement  amplitudes  are  not  needed, 
but  large  acceleration  amplitudes  are  needed.   The  permanent  magnet 
vibration  generator  with  its  large  force  ratings  and  small  moving 
mass  readily  meets  these  requirements.   The  frequency  and  amplitude 
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can  be  varied  simply  by  adjusting  the  oscillator  output.   The  EA  1500 
has  a  frequency  range  from  5  Hz  to  10  kHz  which  is  well  beyond  the 
first  cross  mode  of  operation  and  the  limiting  value  for  this  investi- 
gation. 

The  M-B  Electronics  model  2120MB  power  amplifier  was  selected  to 
complement  the  exciter.   The  distortion  was  less  than  O.77o  at  rated 
output  of  125  watts  in  the  desired  frequency  range  of  operation  and 
less  than  0.3%  at  driving  levels  below  shock.   A  single  unit  provided 
accelerations  to  50g's  for  the  piston  used.   Higher  levels  of  acceleration 
can  be  obtained  by  paralleling  additional  power  amplifiers. 

The  severe  requirements  for  a  very  highly  stable  oscillator  with 
a  precisely  repeatable  frequency  setting  led  to  the  choice  of  the 
General  Radio  Type  1161A-A6C  coherent  decade  frequency  synthesizer  as 
the  primary  frequency  source.   The  synthesizer  has  a  readily  repeatable 
frequency  setting  to  better  than  0.001  Hz.   The  synthesizer  was  checked 
repeatedly  with  a  counter,  and  no  variation  in  the  frequency  output  was 
measured.   The  output  level  control  was  modified  to  permit  finer  adjust- 
ment.  A  General  Radio  Type  1160-P2  Sweep  and  Marker  generator  was  used 
in  conjunction  with  the  synthesizer. 

The  piston  motion  was  monitored  by  the  use  of  a  Endevco  model  2215 
accelerometer  mounted  inside  the  piston.   This  model  will  work  directly 
into  a  VTVM  and  still  provide  satisfactory  frequency  response  up  to 
10  kHz.   It  was  necessary  to  mount  the  accelerometer  on  the  centerline 
of  the  piston  and  directly  to  the  leading  component  of  the  piston  to 
reduce  distortion  and  improve  the  accuracy  of  the  results. 

The  Hewlett-Packard  Electronic   Counter  Model   521C  was   used   to 
monitor    the    frequency    input   and    the    frequencies    from   the  wave   analyzer. 
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No  deviations  were  noted  between  the  set  frequency  and  the  output 
frequency  for  the  synthesizer. 

A  Bruel  and  Kjaer  Type  4136  condenser  microphone  was  mounted 
in  the  rigid  end  so  that  the  diaphragm  was  flush  with  the  end  of 
the  tube.   T  he  small  diameter  (only  \   in.)  was  chosen  to  minimize 
the  disturbance  of  the  rigid  boundary  condition.   The  frequency 
response  of  the  microphone  was  flat  to  within  +  0.5  dB  from  80  Hz 
to  20  kHz,  well  beyond  the  range  of  interest  for  all  five  tubes. 
The  rise  time  of  this  microphone  was  4.5  microseconds. 

All  voltage  measurements  were  taken  with  Hewlett-Packard  400D 
Vacuum  tube  voltmeters,  and  the  waves  shapes  were  observed  using  a 
Tektronic  Type  565  dual-beam  dual-trace  oscilloscope.   The  distortion 
present  in  the  wave  forms  was  analyzed  for  harmonic  content  using 
a  General  Radio  Type  1900-A  wave  analyzer.   This  analyzer  has  a  band 
width  selection  of  3,  10,  and  50  Hz.   Some  Q-curves  were  plotted 
using  data  Equipment  Company  Model  800  AM  X-Y  recorder  and  the 
General  Radio  sweep  and  marker  generator. 

Figure  2  is  a  schematic  of  the  system  used  to  monitor  the 
temperature  of  the  gas.   Because  of  its  linear  resistance  character- 
istic in  the  range  of  interest,  a  Fenwal  Electronics  Inc.  Model 
GA  51P1  thermistor  was  used  and  was  mounted  in  a  small  hole  drilled 
in  the  rigid  end  plate.   The  leads  were  insulated,  and  the  unit  was 
set  in  epoxy  with  the  glass  bead  at  the  surface.   The  epoxy  provided 
an  excellent  seal  and  a  means  of  smoothing  the  end  plate  to  reduce 
distortion.   The  temperature  read-out  was  provided  by  means  of  a 
microammeter ;  when  a  continuous  record  was  needed,  a  strip  recorder 
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was  used.   The  system  was  calibrated  with  the  aid  of  temperature 

baths  of  acetone  and  dry  ice,  ice  and  water,  and  boiling  water. 

With  this  system,  temperature  changes  in  the  tube  of  less  than  0.1  C 

could  be  easily  detected. 

3.2.   Calibration  of  Equipment 

The  experimental  determination  of  P   requires  that  the  absolute 

nn 

calibrations  of  the  microphone  be  known.   In  addition,  the  determination 
of  06  requires  that  the  absolute  calibration  of  the  accelerometer  also 
be  known.   The  accelerometer  sensitivity  is  defined  as  S.  =  V. /A  where 
V.  is  the  peak  amplitude  of  the  accelerometer  output  voltage,  and  the 

microphone  sensitivity  is  defined  as  S_,  =  V^./P  where  V,.,  and  P  are  the 

r  j  M    M         M 

peak  amplitudes  of  the  microphone  output  voltage  and  the  pressure 
respectively . 

The  calibration  of  the  accelerometer  was  accomplished  by  direct 
measurement  of  the  displacement  amplitude  of  the  piston  face  for  a 
range  of  accelerometer  output  voltages.   Displacement  amplitudes  up 
to  0.060  cm  were  measured  with  a  reproducibility  of  0.0002  cm  using 
a  traveling  microscope  and  a  stroboscopic  light.   At  100  Hz,  fifteen 
determinations  of  the  accelerometer  sensitivity  resulted  in 
SA  =  (6-66+  0.03)  x  10"4  volt/(m/sec2). 

The  microphone  was  calibrated  using  a  Bruel  and  Kjaer  Model 
4220  pistonphone  which  produces  a  known  acoustic  pressure  level  of 
124  +  0.2  dB  at  a  frequency  of  250  Hz,   The  sensitivity  was  found 
to  be  SM  =(-76.4  +  0.2)dB  re  volt/^b  =  (1.50  +  0.03)  x  10_J  volt/(Nt/m2) 
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3.3.   Determination  of  the  Attenuation  Constant 

The  necessity  of  having  an  accurately  determined  small  amplitude 
attenuation  coefficient  and  the  reported  wide  variation  in  agreement 
between  the  theoretical  and  the  experimentally  determined  attenuation, 
particularly  in  tubes  containing  standing  waves,  prompted  a  more 
thorough  investigation  of  this  problem. 
3.3.1.   Background 

Beginning  with  Helmholtz  in  1863  (30)^  who  considered  only 
the  effects  due  to  viscosity,  the  causes  of  the  sound  attenuation 
in  cylindrical  tubes  have  been  widely  discussed  and  investigated 
(29,31-42).   Experiments  of  Kundt  led  to  Kirchhoff's  (6)  complete 
theory  for  the  wide  tube  (defined  in  Section  2.2)  which  considered 
both  viscosity  and  heat  conduction.   In  general,  investigators 
report  that  the  Kirchhoff  equation  gives  the  proper  frequency, 
radius,  and  pressure  dependence  but  that  the  experimental  absorption 
is  greater  than  calculated.  Weston  (29)  presented  more  accurate 
expressions  for  the  velocity  and  attenuation  by  including  more  terms 
in  the  expansion  of  the  Bessel  functions.  Weston's  expressions  for 
velocity  and  attenuation  for  the  wide  tube  (which  we  are  considering) 
are: 
Phase  velocity  of  sound: 

cVcfl-'f//;  (2co)^}  (40) 

'•         Hi  — 

-   3,4  34  V  J  ry  f  1  -0> IG  2  -  (t% '  ^ )  ]     cnn  -sec"1 

Attenuation: 

!  •  «  'J 

A'     /  \''s   .-■•)  -  r    •*  -     h"i      "^ !      .'i.,-''  -* 
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where 

y.VT  +  ',i-!)(^'/Ti)'/3---o.l?Ucm/s^  (42) 

The  first  term  of  Eq.  41  is  the  ordinary  wall  absorption,  the  second 
the  ordinary  free  gas  attenuation,  and  the  third  arises  because  the 
sound  energy  is  not  uniformly  distributed  over  the  tube.   The 
numerical  values  here  are  for  dry  air  a  20  C  and  760  mmHg  and  were 
checked  using  constants  obtained  from  National  Bureau  of  Standards 
Circular  564  (43). 

Shields,  Lee  and  Wiley  (44)  present  a  numerical  solution  of  the 
Kirchhoff  equation  and  state  that  Kirchhoff 's  phase  velocity  correction 
is  always  too  large  (giving  values  which  are  too  small),  and  that  the 
calculated  absorption  is  always  too  small.   They  also  give  a  method 
for  determining  the  error  in  the  Helmholtz-Kirchhof f  equations. 

The  work  thus  far  presented  is  for  progressive  waves  in  infinitely 
long  tubes.   It  has  been  stated  by  Shields  (44)  that  neither  the 
Kirchhoff  tube  absorption  and  velocity  equations  nor  the  numerical 
analysis  would  be  expected  to  apply  to  standing  wave  tubes. 

In  considering  resonance  tubes,  there  is  almost  universal  agreement 
on  the  existence  of  excess  experimental  absorption  (44-48)  when  compared 
with  the  absorption  calculated  using  the  equations  for  traveling  waves. 
For  example,  Parker  (45)  reported  an  excess  which  varied  from  6.5% 
for  helium  to  307o  for  argon;  Smith  and  Tempest  (44)  report  30%; 
Henderson  and  Donnelly  (46)  found  107o. 
3.3.2.   Theory  for  Standing  Waves 

For  standing  waves  to  be  generated  in  a  rigid-walled  tube, 
an  integral  number  of  half  wave  lengths  must  be  contained  in  the  tube 
length,  i.e.,  /\    stgb/£m    wnere  m  =  1  >  2 , 3 ,  .  .  .   From  the  fact  that 
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yt    =  c/f 


Jm    =  -TA  (cfei.\  (43) 


where  c  is  the  velocity  of  sound  in  the  fluid  contained  within  the 
tube.  For  a  non-dissipative  medium,  the  infinitesimal  pressure  at 
the  rigid  end  (where  a  =  L)  is  given  by 

?-({&>  A*  on *>*){(&  *<nJ>L)m  (44) 

If  the  medium  is  dissipative,  the  sound  pressure  at  the  rigid 
end  is  still  given  by  Eq.  44,  but  now  k  may  be  considered  as  complex 
in  order  to  account  for  the  dissipation.  We  can  now  write  Eq.  44  as 

which  is  the  same  as  Eq.  28  with  a  90  phase  change  in  the  acceleration. 
And  if  the  absorption  is  small  (i.e.,  0CL<f<'!)>  the  amplitude  of  Eq.  45 
is  given  by  Eq.  30 

p^ippo^/^fjdLf  \   m^hj*  (30) 

This  is  of  the  form  of  a  typical  resonance  curve  if  the  width  of  the 
resonance  curve  is  much  less  than  the  fundamental  frequency.   From  this 
resonance  curve,  the  half-power,  or  the  "down  3dB"  points,  can  be  used 
to  calculate  the  absorption  in  the  resonance  tube  by 

<*e»jB  =  ft/edf*  <46> 

This  method  for  determining  the  attenuation  was  used  by  Parker  and 
by  Henderson  and  Donnelley. 

For  the  resonance  tube,  end  effects  seem  to  be  the  major  additional 
consideration.   If  the  losses  at  the  ends  are  considered  to  be  due  only 
to  heat  conduction  and  are  not  a  function  of  the  viscosity,  then  an 
expression  for  these  losses  may  be  obtained  by  using  the  heat  conduction 
portion  of  Ij    . 
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The  relationship  of  the  viscous  and  heat  conduction  terms  is  given  by 
Eq.  42 

i'-*r  +  a-i)  (*'/*)*  (42) 

where  the  first  term  is  the  viscous  term  and  the  second  term  is  the 
heat  conduction  term.   Assume  that  the  attenuations  are  additive, 
i.e.,  that 

rl  t    (to  fa/)    =  cQ    fQTe  (47 ) 

whereCf.   is  given  by  Eq.  41  and  (X   is  due  to  the  losses  at  the  end. 
t  ;      ~l  e 

Kirchhoff's  expression  for  the  attenuation  due  to  the  walls  alone  is 

K'-a   "    i  ■'  '  ""  OS      (u    ■  c  '  (48) 

which  becomes 

(%cir  -  C.y3(bX<riJ     ,y       /^„  (49) 


when  evaluated  for  dry  air  at  23  C  and  atmospheric  pressure.   Equation 
49  for  tubes  with  a  radius  of  1.125  inches  becomes 

a,(lT-  /.o42 X /o* *jf*    (M"J)>  (50) 

Assuming  that  the  end  has  only  heat  conduction  losses,  then 

ae  i zLCi-omfMu-J' {*><&¥&  '*$= ' &V3^  W*te)      (5i) 

Since  the  heat  conduction  term  is  approximately  one-half  the  viscous 
term  for  dry  air  at  23  C  and  atmospheric  pressure,  the  attenuation  due 
to  the  ends  can  be  related  to  the  attenuation  due  to  the  walls  by 

CCe    ^  &^'ray/3L  (52) 

3.3.3  Experimental 

Attenuation  constants  were  experimentally  measured  for  five 
identical  resonance  tubes  of  varying  length  by  two  methods.   The  first 
consisted  of  measuring  the  width  of  the  infinitesimal-amplitude 


32 


resonance  curve  between  the  ha  If- power  points  and  using  Eq.  46  to 
calculate  the  attenuation  constant.   This  will  be  referred  to  as 
the  Q  method.   The  second  method  consisted  of  measuring  the 
microphone  and  accelerometer  output  voltages  at  the  resonant 
frequency  and  using  Eq.  32  which  can  be  written  in  terms  of 
voltages 

\JH=(fl>  /Ctaxf,  rrr  ft)  ,  £M      SA)   V4  (53) 

A  plot  of  V  vs  V.  ,  for  Q)-&3/»  ,    should  result  in  a  straight  line 

of  slope  (  /O  /(?'  /fft?t)  (S„/S4  )  .   From  this  slope,  QC  can  be  calculated 

f° ' '   -tYf  MA  exp 

if  S^,  and  S.  are  known.   This  method  will  be  referred  to  as  the  A/P 
M      A 

method. 

3.3.4.   Comparison  of  results 

Attenuation  constants  were  calculated  for  all  five  tubes  using 

Eq.  49  and  Shield's  corrections  for  dry  air  at  23  C  and  760  mmHg. 

These  values  for  CX(l,j  axe.   compared  with  the  experimentally  determined 

6"> ''.'  for  the  fundamental  modes  in  Table  1.   The  experimental  results 

for  all  five  tubes  by  both  methods  are  compared  with  the  theory  in  the 

range  of  interest  in  Figs.  3  and  4  where  the  log  of  the  attenuation  is 

plotted  vs  the  log  of  the  frequency.   There  is  almost  perfect  agreement 

between  the  calculated  and  the  experimentally  observed  attenuation 

constants  by  the  Q  method.   This  is  true  through  the  50th  harmonic 

for  the  9'  tube  and  the  40th  harmonic  for  the  6'  tube.   At  these 

points  side  modes  were  calculated  to  be  present  and  were  observed  by 

distortion  in  the  Q-curves. 

Except  for  the  7  in.  tube,  the  effect  of  the  end- loss  correction 

QC     is  too  small  to  see  on  the  graphs.  As  a  check  on  the  end- loss 

approximation,  if  uC  =  <X   ,  +GC    and  GC  in  Table  1  are  compared, 

x      cal    c        exp 

it  is  seen  that  the  assumption  that  the  end  losses  are  due  solely 
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to  heat  conduction  is  slightly  low.   If  the  same  losses  are 

assumed  for  the  ends  as  for  the  side  wall,  then  an  R/L  relationship 

would  hold.   Again  compare  (X  m  and  r£  in  Table  1,  and  it  can 

T       exp 

be  seen  that  this  correction  is  too  large  indicating  that  the 
end  correction  is  somewhere  between  these  two  extremes.   If  an 
empirical  relationship  of  R/2L  is  taken,  then  the  agreement  is 
better  than  1%   for  all  but  the  seven  inch  tube. 

The  results  of  these  investigations  indicated  that  the  Kirchhoff 
equation  (using  Shield's  corrections  and  properly  corrected  for  end 
losses)  can  be  used  to  determine  the  attenuation  in  standing  wave 
tubes  of  high  Q  where  the  boundaries  approach  ideal  boundaries  and 
where  the  length  of  the  tube  is  much  greater  than  the  radius  (on 
the  order  of  ten  to  one) . 

The  results  further  indicate  that  both  the  Q  method  and  the  A/P 
method  are  valid  methods  for  determining  the  attenuation.   The  Q 
method  has  the  advantage  of  requiring  very  little  data  and  data 
processing  to  obtain  the  attenuation  when  a  highly  accurate  frequency 
source  is  available.  However,  the  acceleration  had  to  be  closely 
controlled,  or  an  attenuation  constant  greater  than  actually  existed 
was  found.   This  change  in  the  impedance  of  the  driver  near  resonance 
could  account  for  some  of  the  problems  in  the  past.   The  A/P  method 
is  reliable  and  useful  when  the  frequency  is  not  critical  though  the 
accelerometer  and  microphone  must  be  accurately  calibrated. 

The  effects  of  temperature  cannot  be  minimized.   A  3  C  change 
makes  approximately  a  l"/0   change  in  Kirchhoff  s  constant  for  dry  air 
as  given  by  Eq.  42  and  a  resulting  change  in  (X     .   Also  of  interest 
is  that  a  0.6  C  temperature  change  will  result  in  a  0.1  Hz  change 
in  U\,    for  the  six  foot  tube.   Therefore,  it  was  important  that  the 
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temperature  of  the  system  be  closely  monitored,  particularly  when 

the  resonance  tube  was  being  excited  at  very  high  levels. 

3.4.   Determination  of  the  Resonance  Frequency 

The  resonance  frequency  was  determined  by  finding  the  frequency 

at  which  P.. ,  was  a  maximum  for  a  fixed  A  .   Since  £&■*>   was  defined  on 
11  o  ' 

the  basis  of  the  first-order  solution,  A  was  set  low  enough  that 
finite-amplitude  effects  were  negligible.   Repeated  measurements 
indicated  that  -J-*    could  be  determined  to  better  than  0.01  Hz 
providing  consideration  was  given  to  temperature.  As  a  check  on  (£2^    , 
the  frequencies  of  the  "down  3  dB1'  points  (which  could  be  determined 
to  within  0.001  Hz)  were  determined,  and  the  mid- frequency  was  compared 
with  00 .      The  agreement  was  within  0.005  Hz. 

To  arrive  at  a  value  of  eOp    characteristic  of  the  conditions 
actually  prevailing  at  the  time  that  the  finite-amplitude  data  were 
collected,  the  experimental  procedure  was  designed  so  that  W^  and  the 
temperature  were  determined  before  the  run  and  after  the  run  and  the 
temperature  was  recorded  throughout  the  run.   This  technique  insured 
that  OJ^  was  known  to  within  approximately  0.02  Hz  for  all  runs. 
3.5   System  Alignment 

It  was  necessary  that  the  piston  be  very  carefully  aligned  before 
each  run  or  distortion  was  introduced  into  the  input  waveform.   The 
system  was  designed  so  that  the  resonance  tube  was  fixed  and  the  shaker 
was  mounted  on  an  vibrational  insulated  base  plate  that  was  provided 
with  three  dimensional  adjustment. 

The  alignment  procedure  consisted  essentially  of  adjusting  the 
piston  in  the.  tube  until  it  was  free  to  move.   The  piston  was  then 
excited  at  a  low  level,  and  the  accelerometer  output  was  monitored  on 
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one  trace  of  a  dual  trace  oscilloscope.   A  sine  wave  from  the 
frequency  synthesizer  was  displayed  on  the  other  trace.   Fine 
adjustment  was  made  to  the  shaker  base  until  the  two  wave  forms 
appeared  to  be  identical.   The  excitation  level  was  then  increased 
to  the  maximum  level  desired  and  additional  fine  adjustment  was 
made  as  necessary. 

The  distortion  present  was  determined  by  measuring  with  a 
wave  analyzer  the  percentage  of  harmonic  content  relative  to  the 
fundamental  in  the  accelerometer  output.   Due  to  the  impedance 
mismatch  between  the  accelerometer  and  the  wave  analyzer,  it  was 
necessary  to  use  some  form  of  impedance  matching,  and  a  Burr-Brown 
Model  100  preamplifier  was  used  because  it  had  a  low  noise  and 
provided  excellent  isolation  characteristics. 

By  using  the  visual  method  of  waveform  comparison  for  alignment, 
it  was  found  that  the  second  harmonic  distortion  of  the  acceleration 
waveform  could  be  kept  below  0.27o  of  the  fundamental  and  the  third 
harmonic  distortion  could  be  kept  below  0.37o  of  the  fundamental. 
All  other  harmonics  were  masked  by  the  noise  level. 
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4.   Measurements  and  Analysis 

Three  measurement  procedures  were  employed  to  compare  the 
experimental  results  with  the  theory.   These  procedures  provided 
information  to  compare  the  pressure  waveforms,  Q-curves,  and  a. 
Fourier  analysis  of  the  waveforms  with  these  quantities  predicted 
by  theory. 

4.1.  Pressure  Waveform 

The  pressure  waveforms  for  various  sound  pressure  levels  were 
compared  with  those  predicted  by  the  computer  program  FINAMPI.   The 
program  FINAMPI  computes  and  graphs  the  predicted  pressure  waveform 
at  the  rigid  end  using  all  terms  in  the  perturbation  solution  through 
the  sixth  order  with  input  variables  of  , u) '/.(."; and  M. 

The  measurement  procedure  consisted  of  first  setting  A  and  the 
frequency  at  either  &)«>    or  at  some  desired  frequency  off  resonance. 
The  pressure  waveform  was  then  monitored  on  one  channel  of  a  dual- 
beam,  dual-trace  oscilloscope.   A  sine  wave  from  the  frequency 
synthesizer  was  displayed  on  the  other  channel  for  reference  and 
distortion  determination.   Photographs  were  made  of  the  oscilloscope 
trace  for  comparison  purposes. 

The  analysis  consisted  simply  of  visually  comparing  the  photographed 
pressure  waveforms  at  the  various  SPL's  and  frequencies  with  those 
predicted  by  the  computer  program. 

4.2.  Fourier  Analysis 

A  Fourier  Analysis  of  the  pressure  waveform  was  made.   The  same 
procedure  was  used  here  as  with  the  waveform  analysis.   In  addition  the 
harmonic  content  was  measured  as  for  the  Q-curves. 
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The  analysis  consisted  of  applying  a  grid  to  the  photograph 
and  taking  the  amplitude  of  64  equally  spaced  points  in  one  cycle. 
These  data  points  were  used  as  the  input  for  the  computer  program 
FOUANAL . FOUANAL  takes  these  data  points  and  does  a  Fourier  snythesis 
on  the  wave  shape.   The  amplitudes  of  the  Fourier  components  were 
then  compared  with  the  measured  values  of  the  harmonics. 
4.3.   Q-Curves 

It  is  possible  to  relate  the  response  of  the  resonance  tube  to 
the  frequency  response  of  an  electrical  network.   In  this  case  the 
pressure  curves  of  the  resonance  tube  would  be  analogous  to  the  "Q" 
response  curves  of  the  electrical  network.   Ruff  has  developed  an 
expression  for  the  amplitude  of  each  harmonic  of  the  pressure  waveform 
as  a  function  of  ZiiU)  and  M.   The  procedures  outlined  here  provided 
the  means  for  displaying  the  harmonic  content  of  the  pressure  wave 
as  a  function  of  frequency  near  resonance  for  a  desired  value  of  A  . 

The  measurement  procedure  consisted  of  first  setting  the  value 

of  A  and  then  measuring  with  a  wave  analyzer  the  pressure  of  the  nth 

harmonic,  P   ,  at  the  rigid  end  as  a  function  of  frequency.   The 
'   nn  °  i  j 

frequency  was  swept  from  well  below  the  half-power  point  to  well  above 
the  half- power  point  for  the  fundamental  in  steps  of  0.02  Hz.   For  the 
other  harmonics  through  the  sixth,  the  frequency  was  swept  from  just 
above  the  noise  level  through  the  area  of  interest  back  to  just  above 
the  noise  level  in  steps  of  0.01  Hz. 

The  data  were  analyzed  by  first  plotting  the  various  P   in  te: 

J  J  c  °  nn 

of  voltage  versus  the  frequency.   The  data  were  then  normalized  using 

the  magnitude  of  the  leading  term  of  the  particular  harmonic  taken 

from  the  computer  program  QCURVE  and  QCURC.   This  was  done  to  facilitate 


38 


comparison  with  the  normalized  graphs  that  were  drawn  by  the  computer. 
The  fundamental  resonance  frequency  was  determined  and  used  as  the 
reference  frequency  for  comparison  with  the  computer  graphs.   The 
normalized  data  points  for  each  harmonic  were  plotted  and  compared 

to  the  computer  graphs  provided  by  the  two  programs,  QCURVE  and  QCURC. 

QCURVES  is  a  program  that  uses  the  perturbation  solution  and  is 
designed  to  compute  and  graph  the  resonance  response  of  each  harmonic 
using  terms  of  the  sixth  order  as  the  fundamental  is  swept  through 
resonance.   QCURC  is  a  program  similar  to  QCURVES  except  that  it  uses 
the  leading  terms  of  each  harmonic  in  the  perturbation  solution. 

The  major  points  considered  were  the  amplitude  of  the  harmonic, 
the  frequency  at  which  the  maximum  amplitude  for  each  harmonic  occurred, 
and  the  bandwidth  of  the  various  harmonics. 
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5.   Results 

In  this  section  the  experimental  results  will  be  compared  with 
the  theoretical  results  obtained  by  Ruff. 
5.1.   Pressure  Waveforms 

The  first  comparisons  made  will  be  of  the  pressure  waveforms 
(Figs.  5-8).   In  these  figures  the  large  drawing  is  of  the  waveform 
predicted  by  the  perturbation  approach  (program  FINAMPI)  and  the 
smaller  inset  is  a  drawing  from  an  oscilloscope  photograph  of  the 
actual  pressure  waveform.   The  dashed  line  in  both  drawings  is  a 
sine  wave  which  is  included  as  a  reference. 

In  the  low  Mach  number  region  which  includes  Figs.  5  and  6,  it 
appears  that  the  theory  and  experiment  are  in  very  close  agreement. 
In  this  region  there  is  relatively  little  finite-amplitude  distortion, 
and  good  agreement  is  expected.  When  the  Mach  number  exceeds  a  value 
of  about  M  =  0.005,  the  finite-amplitude  distortion  reaches  a 
significant  level.   However,  Fig.  7,  which  is  for  a  Mach  number  of 
M  =  0.007  indicates  that  the  theoretical  and  experimental  results 
are  still  in  surprising   agreement.  As  the  Mach  number  is  increased 
into  the  realm  between  significant  finite-amplitude  effects  and  the 
region  of  shock,  the  theory  begins  to  breakdown  as  is  evidenced  in 
Fig.  8.   This  figure  is  for  M  =  0.009,  and  it  is  seen  that  the 
theoretical  waveform  displays  details  which  are  not  observable  in  the 
experimental  waveform. 

It  is  evident  that  the  theory  becomes  invalid  when  it  is  extended 
into  the  region  of  relatively  high  Mach  number. 
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5.2.  Fourier  Analysis 

In  addition  to  the  foregoing,  the  experimentally  determined 
values  of  the  amplitude  of  each  harmonic  may  be  compared  with  the 
results  from  the  program  FOUANAL. 

The  results  of  the  comparison  with  the  program  FOUANAL  are 
included  in  Table  2.   FOUANAL  uses  the  data  taken  point  by  point 
from  an  oscilloscope  photograph  and  computes  the  Fourier  coefficients 
of  the  various  terms  of  a  Fourier  expansion  of  the  waveform.   The 
small  size  of  the  photographs  and  the  thickness  of  the  trace  line 
in  the  photographs  were  such  that  only  qualitative  results  were 
expected. 

FOUANAL  also  computes  the  phase  angle  associated  with  each 
harmonic.   These  results  are  included  in  Table  3. 

5.3.  Q-Curves 

Q-curves  for  three  Mach  numbers  (M  =  0.004,  M  =  0.005,  M  =  0.009) 
were  experimentally  obtained  and  graphed  for  comparison  with  the  Q- 
curve  program  results  for  the  same  Mach  numbers.   These  computer  program 
results  indicated  by  a  dashed  line  for  QCURVES  and  a  solid  line  for 
QCURC  are  contained  in  Figs.  9  through  26.   Several  things  should  be 
mentioned  concerning  these  graphs.   First,  the  three  Mach  numbers  chosen 
cover  the  range  from  very  slight  distortion  to  significant  distortion 
and  should  indicate  the  usable  extent  of  the  theory.   Second,  in  every 
figure,  the  Q-curve  results  from  QCURVES,  lies  below  the  results 
obtained  using  QCURC  (the  first  term  method).   The  data  points  in 
general  fall  very  close  to  the  results  of  QCURC  except  for  very  high 
Mach  numbers  where  the  predicted  amplitude  of  the  higher  harmonics  is 
much  larger  than  the  experimentally  determined.   Another  interesting 
fact  is  that  the  same  harmonics  for  each  Mach  number  have  exactly  the 
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same  Q-curve  using  the  first  term  method.   This  is  not  surprising 
when  it  is  remembered  that  in  the  normalization  scheme  the  dependence 
of  the  leading  terms  will  not  change  for  different  Mach  numbers. 
5.4.   General 

It  would  have  been  interesting  to  obtain  the  same  results  for 
the  other  tubes  used  in  the  attenuation  measurements,  but  time  did 
not  permit.   However,  some  observations  that  were  made  are  worthy  of 
note.   The  distortion  of  the  pressure  waveform  developed  in  the  same 
way  for  all  five  tubes;  however,  in  the  formation  of  shock  slight 
differences  were  noted  and  were  more  prominent  in  the  short  tube. 

A  sequence  of  photographs  of  the  pressure  waveforms  for  the 
14  in.,  6  ft.,  and  9  ft.  tubes  at  sound  pressure  levels  of  155,  157, 
and  160  dB  were  taken  and  are  shown  in  Figs.  27,  28,  and  29. 

The  wave  shapes  are  very  similar;  however,  there  are  certain 
noticable  differences.   The  small  spur  at  the  base  of  the  shock  is 
present  in  all  waveforms,  but  is  more  prominent  in  those  for  the  6  ft. 
tube.   This  point  needs  further  investigation  to  determine  whether 
the  pressure  actually  decreases  immediately  before  the  shock  or  if 
the  spur  is  due  to  some  property  of  the  instrumentation.   The  steepness 
of  the  shock  front  increases  as  the  resonant  frequency  increases.   The 
grass  in  the  pressure  waveform  increases  with  an  increase  in  resonant 
frequency  and  the  drive  level.   This  is  a  result  of  the  "ringing"  in 
the  piston  introduced  by  the  Shockwave  banging  against  the  piston  face. 
The  "grass"  is  also  more  prominently  displayed  at  the  higher  frequencies 
due  to  the  relative  time  scales  and  the  higher  drive  levels  required  to 
attain  the  same  sound  pressure  level.  At  very  high  drive  levels,  the 
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shock  wave  striking  the  piston  can  be  observed  in  the.  accelerometer 
output  from  the  piston,  and  its  position  can  be  shifted  by  small 
changes  in  the  drive  frequency. 
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6.   Conclusions 

Experimental  investigation  has  shown  that  the  theory  accurately 
predicts  the  position  relative  to  the  fundamental  at  which  the  maximum 
amplitude  of  each  harmonic  occurs  throughout  the  range  from  very  low 
to  very  high  Mach  numbers.   The  investigation  has  also  shown  that  in 
the  region  of  relatively  low  Mach  number  the  theory  accurately  predicts 
the  amplitude  of  the  first  six  harmonics  that  comprise  the  pressure 
waveform.   The  fact  that  the  third  harmonic  is  slightly  high  at  the 
lower  Mach  numbers  can  be  accounted  for  in  part  by  the  third  harmonic 
distortion  in  the  input.   At  high  Mach  numbers  the  predicted  amplitudes 
of  the  higher  harmonics  are  much  greater  than  the  experimentally 
observed  ones.   At  this  writing  there  is  no  definite  explanation  of 
this . 

The  reconstructed  pressure  waveforms  drawn  by  the  computer  are  in 
very  close  agreement  with  the  experimentally  observed  ones  at  relatively 
low  Mach  numbers.   When  the  finite-amplitude  effects  become  pronounced, 
the  theory  does  not  properly  predict  the  spur  that  forms  just  prior  to 
the  shock  formation;  however,  the  theory  does  indicate  a  definite 
distortion  of  the  waveform  in  this  region. 

The  Fourier  analysis  of  the  pressure  waveform  shows  good  agreement 
with  the  experimentally  determined  amplitudes  at  all  Mach  numbers. 
The  experimental  phase  angles  show  an  excessive  scatter,  and  the 
agreement  with  the  theoretical  values  is'  only  qualitative.   This  appears 
to  be  the  only  place  where  the  experimental  results  are  less  reliable 
than  the  theory.   Since  the  major  part  of  the  distortion  takes  place 
in  a  relatively  small  portion  of  the  waveform,  the  accuracy  of  the 
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Fourier  analysis  could  be  improved  through  the  use  of  more  data  points 
in  the  region  of  the  developing  shock.   It  is  suggested  that  further 
work  along  this  line  be  postponed  until  the  oscilloscope  waveform  can 
be  more  closely  linked  with  the  computer  input. 

In  addition,  the  dependence  of  the  theory  on  the  infinitesimal- 
amplitude  attenuation  coefficient  has  been  demonstrated. 

In  general,  then  the  statement  that  the  perturbation  theory  can 
adequately  describe  the  effects  of  finite-amplitude  standing  waves 
in  rigid-walled  tubes  has  been  shown  to  be  valid  provided  the  Mach 
number  is  such  that  the  distortion  is  not  too  great.   This  of  course 
corresponds  to  the  pre- shock  area,  and  even  then  precludes  much  of 
that  area  in  which  significant  distortion  is  evident. 
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ATTENUATION   CONSTANTS   AND  END-EFFECT  CALCULATIONS   FOR 
FIVE    RIGID-WALLED   TUBES   CONTAINING   STANDING  WAVES 

TABLE    1 
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